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Nanoconfined catalytic Angstrom-size motors 
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Canada 

Self-propelled chemically powered synthetic micron and nano-scale motors are being intensively studied be¬ 
cause of the wide range of potential applications that exploit their directed motion. This paper considers 
even smaller Angstrom-size synthetic motors. Such very small motors in bulk solution display effects arising 
from their self-propulsion. Recent experiments have shown that small-molecule catalysts and single enzyme 
molecules exhibit properties that have been attributed to their chemical activity. Molecular dynamics is used 
to investigate the properties of very small Angstrom-size synthetic chemically powered sphere-dimer motors 
in a simple atomic-like solvent confined between walls separated by distances of tens of nanometers. Evi¬ 
dence for strong structural ordering of the motors between the walls, which reflects the finite size of solvent 
molecules and depends on solvent depletion forces, is provided. Dynamical properties, such as average motor 
velocity, orientational relaxation, and mean square displacement, are anisotropic and depend on the distance 
from the walls. This research provides information needed for potential applications that use molecular-scale 
motors in the complex confined geometries encountered in biology and the laboratory. 


I. INTRODUCTION 

Synthetic self-propelled motors that convert chemical 
energy from their environment into directed motion are 
examples of active objects with distinctive and useful 
properties. Much of the interest in such motors stems 
from potential applications that exploit their directed 
motion to carry out functions involving cargo transport, 
analogous to the active transport tasks performed by 
molecular motors in the cell . 1-5 

In many circumstances, motors operate in complex en¬ 
vironments or systems confined by boundaries. Examples 
include motors used in microfluidic devices 6 and poten¬ 
tial in vivo applications such as targeted drug delivery 7,8 . 
Also, in many experiments currently being carried out, 
motors reside near a surface as a result of gravitational 
or other forces, and the interactions of motors with sur¬ 
faces can lead to interesting dynamical effects . 9-12 For 
these reasons, it is important to assess the influence of 
boundaries and confinement on motor dynamics. 

The effects of confinement on motor motion have been 
investigated previously. The velocity of a spherical self- 
propelled particle with a point-like catalytic site con¬ 
fined by a spherical wall has been computed analyt¬ 
ically . 13 Micron-scale Janus motors operating by self- 
diffusiophoresis have been studied experimentally and 
numerically in circumstances where they collide with a 
wall, are confined to a circular pore, and move through a 
triangular lattice of obstacles forming a patterned envi¬ 
ronment . 14 The dynamics of a self-diffusiophoretic Janus 
particle close to a hard wall has been shown to ex¬ 
hibit reflection, steady sliding, and hovering . 15 A study 
of such Janus particles near a wall has been carried 
out using continuum hydrodynamics . 16 The rectification 
of a Brownian Janus particle in a triangular channel 1 ' 


a ) Electronic mail: pcolberg@chem.utoronto.ca 
^Electronic mail: rkapral@chem.utoronto.ca 


and the escape of an ellipsoidal Brownian Janus par¬ 
ticle from a two-dimensional sinusoidally corrugated or 
square-shaped pore have been simulated . 18 Studies of 
confined self-propelled rods have shown that they tend to 
reside near surfaces . 19 The accumulation of self-propelled 
Brownian spheres near a hard wall has been studied 
analytically and simulated using multi-particle collision 
dynamics . 20 Similarly, ensembles of self-propelled col¬ 
loidal rods aggregate and structurally order on confining 
walls . 21 

In a related but somewhat different context, the effects 
of hydrodynamic interactions on the confined motions of 
swimmers and biological organisms have also been in¬ 
vestigated. Since biological organisms function in envi¬ 
ronments containing obstacles and boundaries of vari¬ 
ous types, such studies are essential in order to obtain 
a full understanding of swimming motions in a biolog¬ 
ical context. Studies of hydrodynamic effects on self¬ 
propulsion of both synthetic swimmers and biological or¬ 
ganisms near boundaries and in complex media have been 
carried out . 21-26 Boundaries can influence the collective 
behavior of swimmers and the collective motions of swim¬ 
mers confined to the gap between two planar surfaces 
have been observed . 24,25 

Several factors must be considered in order to un¬ 
derstand the effects of confinement on chemically pow¬ 
ered motors. These include direct interactions with the 
boundaries, motor geometry, the influence of boundaries 
on hydrodynamic flows, and the nature of chemical gra¬ 
dients near walls. At the micron to nanometer scales, 
self-propulsion is influenced by thermal fluctuations. Ori¬ 
entational Brownian motion limits the duration of the 
ballistic regime where directed motion can be observed 
and leads to diffusive dynamics on long time scales, al¬ 
beit characterized by diffusion coefficients that reflect the 
underlying active motion. Interactions with boundaries 
can influence orientational Brownian motion and modify 
motor behavior. 

In this paper, we consider a regime of motor dynamics 
and confinement that differs from that in earlier inves- 
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tigations. We focus on very small Angstrom-scale mo¬ 
tors confined in thin nanometer-scale layers between two 
walls. This study is prompted by several observations. 
Recent experiments have demonstrated that catalytically 
active 5 A organometallic molecules 27 as well as single en¬ 
zyme molecules 28,29 display enhanced diffusion compared 
to the same systems in the absence of chemical activity. 
In addition, experimental studies of very small (30 nm) 
synthetic self-propelled Janus particles have also found 
enhanced diffusion coefficients . 30 Although the precise 
mechanism responsible for diffusion enhancement is still 
a matter of debate for the molecular systems, the ob¬ 
served phenomena are correlated with chemical activ¬ 
ity. Such enzymatically induced motion has been re¬ 
viewed recently . 31 In addition to these experimental re¬ 
sults, molecular dynamics simulations of chemically pro¬ 
pelled Angstrom-size sphere-dimer motors that operate 
by self-diffusiophoresis have shown that these motors can 
move distances of one to several times their length before 
they reorient and have enhanced diffusion coefficients . 32 

For such small-scale motors, whose size is compara¬ 
ble to the molecules comprising the environment, solvent 
structural effects play an important role and fluctuations 
dominate their behavior. Nevertheless, the signatures 
of self-propulsion are clearly evident in observable prop¬ 
erties such as diffusion coefficients. Enzymes carry out 
their functions in the complex cellular environment and, 
potentially, very small synthetic motors will find appli¬ 
cations on the cellular level and other confined systems. 
It is then important to consider effects due to the con¬ 
finement on these very small motors. 

System-specific features of both the motor and solvent 
molecules can determine the structural and dynamical 
properties of motors with small molecular-scale dimen¬ 
sions. Some general properties are common to all such 
motors. Small motors reorient very rapidly, for example, 
nanometer-scale motors have reorientation times of the 
order of nanoseconds or picoseconds, compared to times 
of the order of seconds for micron-sized motors. This time 
scale disparity will have an effect on the relative impor¬ 
tance of ballistic and diffusive motions. For small mo¬ 
tors, solvent structural correlations can influence effec¬ 
tive motor-wall forces in confined systems, as well as the 
concentrations of fuel species in the motor vicinity. Al¬ 
though the quantitative form these effects take depends 
on the system under study, their presence is a common 
feature of molecular-scale motors. In this study, we con¬ 
sider a simple dimer motor made from catalytic and non- 
catalytic spheres in a dense atomic-like solvent in order 
to investigate some of the general aspects of solvent and 
confinement effects on very small motors. 

The outline of the paper is as follows: In Sec. II, we 
describe the model for the motor and solvent molecules, 
along with their interactions with the confining walls. 
Section III presents results on static structural proper¬ 
ties, in particular, the nature of the motor position and 
its angular distributions perpendicular and parallel to the 
walls. Dynamical properties, including dimer velocity, 


orientational correlations, and mean square displacement 
(MSD), are discussed in Sec. IV. The conclusions of the 
study are given in Sec. V. 

II. CONFINED SELF-PROPELLED MOTORS 

We consider motors that achieve propulsion by uti¬ 
lizing catalytic reactions through a self-diffusiophoretic 
mechanism . 5,33-35 In this mechanism, a chemical reac¬ 
tion on a portion of the motor leads to asymmetric dis¬ 
tributions of reactants and products in the motor vicin¬ 
ity. If the reactants and products interact with the motor 
through different intermolecular potentials, this gives rise 
to a body force on the motor. Since no external forces 
are applied to the system, the entire system, the motor 
and its environment, are force-free, and momentum con¬ 
servation leads to fluid flows in the environment. This 
mechanism is responsible for the directed self-propelled 
motion of such motors. 

The most widely studied motors of this type are Janus 
colloids where one face is catalytic and the other chemi¬ 
cally neutral. Suppose that the chemical reaction A —► B 
occurs on the catalytic face and the A and B molecules 
interact with the Janus particle through t/j ,a and C/j.b 
intermolecular potentials. Using the equations of contin¬ 
uum hydrodynamics, supplemented with boundary con¬ 
ditions that account for reactions and fluid velocities on 
the surface of the Janus particle, the velocity of the motor 
is given by 

V = ( k B T/r,)Aj(z • V,c B (rj,0)) 5 , (1) 

where z is a unit vector from the noncatalytic to catalytic 
faces of the Janus particle, CB(rj, 6) is the concentration 
of species B on the surface as a function of the angle 0, 
(• • • )s denotes the average over the spherical surface, and 

nOO 

Aj = / dr r[e-^ J,B W _ e -/»tfj.A(r)], 

Jo 

with /? = 1 /(/cbT). From this expression, we see that the 
velocity depends on the surface average of the gradient 
of the concentration along the surface, the difference in 
the interaction potentials as determined by Aj and 77 , the 
viscosity of the solution. These are features of all motors 
that operate by diffusiophoretic mechanisms. 

Motor geometry is a factor to consider when dealing 
with wall interactions, and instead of spherical Janus mo¬ 
tors, we investigate the dynamics of sphere-dimer motors 
that are simple examples of more complex, non-spherical 
motors. A sphere-dimer consists of linked catalytic (C) 
and noncatalytic (N) spheres. In this geometry, the cat¬ 
alytic region is restricted to one sphere . 36 The directed 
motion of the motor is again due to a diffusiophoretic 
mechanism. Experimental 11 and simulation 37-39 studies 
of the motions of sphere-dimer motors in bulk solution 
have been performed. The analytic continuum theory 
for the velocity of a sphere-dimer motor is considerably 


3 


more involved due to the motor geometry, but it may be 
obtained using a bispherical coordinate system . 40,41 The 
theory predicts that the velocity of the motor along the 
dimer bond is given by a form analogous to that for the 
Janus particle, 


V = ( k B T/ri)AsDK , (3) 

with K given by a complicated but explicit expression 41 
that plays the same role as the surface average of the 
concentration gradient for the Janus particle. If the in- 
termolecular potentials only differ for interactions of A 
and B with the noncatalytic sphere, A sd has the form 

/•OO 

A sd = dr r[e-^ N,B ( r) - e -^ N,A ( r )]. (4) 

Jo 

Since we are interested in very small Angstrom-scale 
motors, continuum theory is not able to describe all rele¬ 
vant effects. As noted earlier, fluctuations are especially 
important and the finite size of the solvent molecules is 
comparable to that of the motor, both of which violate 
the assumptions of a continuum description of the sol¬ 
vent. Instead we adopt a full molecular dynamics de¬ 
scription of a self-propelled Angstrom-scale motor con¬ 
fined between two parallel walls. 

In more detail, the sphere-dimer motor we consider 
consists of a catalytic C sphere and a noncatalytic 
N sphere, with diameters ac = 2a and ctn = 4a, re¬ 
spectively, linked by a rigid bond of length d. While 
various dimer-sphere sizes can be chosen, and the results 
will depend on these values, we have selected a sphere- 
diameter ratio of (tc/on = 0.5 that has been shown to 
yield the highest propulsion velocity both experimen¬ 
tally 11 and theoretically 41 . The solvent consists of N s 
structureless A and B particles with diameter a and mass 
m. The solvent number density is g s = N s /V = 0.8a -3 , 
which corresponds to a dense fluid. The sphere dimer 
masses, M m = |p s to(ctq + a ^), are chosen to make the 
dimer neutrally buoyant. All particles interact through 
a shifted, truncated Lennard-Jones potential, Vij(r) = 
eij{4[(cr i j/r ) 12 - (< 7ij/r ) 6 ] + 1 } for r < \/2aij and zero 
otherwise. Here r is the distance between the centres 
of a pair of particles, ctca = acB = \ (Ac + cr) and 
cjna = unb = \ An + a) for pairs of dimer sphere and 
solvent particle, and ctaa = ctab = cbb = cr for pairs 
of solvent particles. The interaction energy is e for all 
pairs apart from NB pairs, where cnb = O.le or lOe. The 
temperature of the system is k^T/e = 1 . 

The system is contained in a slab with edge lengths 
L x = L y = 50u and L z = 10,15,..., 50cj with 
periodic boundary conditions in X and Y directions 
and two planar walls at Z = 0 and L z . The dimer 
spheres and solvent particles interact with the walls 
via a repulsive 9-3 Lennard-Jones potential, VAfC) = 
e w s{(3\/3/2) [(a w s/C) 9 -(a wS /C) 3 ]+1} for < < \/3a w s, 
where ( is the distance to the closest wall, and zero other¬ 
wise. The wall interaction parameters are a w a = u w b = 



FIG. 1. Cross section of sphere-dimer motor with spheres C 
(blue) and N (vermilion) immersed in A (sky blue) and B 
(orange) solvent particles confined between two planar walls 
a distance Lz = 20 apart. 


cr ; a w c — (ac T a)/2, a W N — (un T a)/2, c w a — CwB — e, 
and e w c = c w n = e or 2 e. 

The C sphere catalyses the reaction A —> B, where A is 
converted to B with unit probability when it lies within 
a distance \/2 acB from the C sphere. The system is 
maintained in a nonequilibrium steady state through a 
bulk back-reaction B —> A with rate 10~ 3 r _1 outside 
of the C- and N-sphere interaction zones. The dimer 
bond length is taken to be d = \/2 (ucb + c^nb), which is 
chosen to ensure energy conservation in the presence of 
reactions. 

Full molecular dynamics 42 is used to follow the motions 
of the motor and solvent and capture their structural and 
dynamical properties . 43 The model neglects the structure 
of the solvent molecules, but nevertheless it is suited to 
reproduce structural and dynamical effects occurring at 
the Angstrom scale. As discussed in the Introduction, 
applications to experimental systems should take into ac¬ 
count the detailed structural properties of the motor and 
solvent under investigation. The velocity-Verlet integra¬ 
tion time step for the molecular dynamics simulations is 
10~ 3 r. Simulation results are reported in dimensionless 
units with distance given in units of a, mass in units of 
m, energy in units of e, and time in units of r = a\/me -1 . 

Parameters were chosen to model a dense fluid Argon- 
like solvent. Given Argon 44 values of a = 0.34 nm, 
e = 120Kfc B , and m = 39.95 u and r = 2.15 ps, we can 
assign physical values to our Angstrom-scale motor sim¬ 
ulations. In these units the sphere-dimer monomers have 
radii of 0.34 nm and 0.68 nm for the C and N spheres, re¬ 
spectively. The end-to-end length of the sphere-dimer 
motor is 2.55 nm. The separation between the walls 
varies from 3.4 nm to 17nm. 

The continuum theory for sphere-dimer motors that 
operate by self-diffusiophoresis has been carried out 40,41 
and the fluid flow fields associated with propulsion have 
been determined 41 . Motors with cnb = 0.1 and small 
dimer bond lengths move in the direction of the C sphere 
and the far-held fluid how is characteristic of a puller. 
For cnb = 10, the motor is propelled in the direction of 
the N sphere and the far-held fluid how is characteristic 
of a pusher. Both types of motor were investigated in 
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confinement between two planar walls, and in bulk with 
periodic boundary conditions in all directions. 

To isolate effects that are due to the self-propulsion 
of the motor, results for an inactive dimer (cnb = 1 ) 
are also presented. Inactive dimers were constructed by 
turning off the chemical reactions that power the mo¬ 
tor and maintain the system in a nonequilibrium steady 
state. The effect of the wall potential on the dimer is il¬ 
lustrated by varying the repulsiveness between walls and 
spheres, denoted in the following as the weaker (e w g = 1) 
and the stronger wall potential (e w s = 2). 

An instantaneous configuration of the sphere dimer 
and surrounding A and B solvent species confined be¬ 
tween two walls is shown in Fig. 1. In this figure one 
can see the relative sizes of dimer monomers and solvent 
particles, as well as the inhomogeneous distribution of 
reactants and products near the dimer and wall, which 
will play an important role in the subsequent discussion. 


III. STRUCTURAL ORDERING OF DIMER MOTOR IN 
CONFINED GEOMETRY 

We consider the structural ordering of the sphere- 
dimer motor when it is confined between the two walls. 
Sample trajectories of an inactive dimer, and puller and 
pusher motors are shown in Fig. 2. The dimer is ini¬ 
tially placed along the middle plane between the walls. 
Depending on the absence or presence of propulsion, the 
dimer takes a long or short time, respectively, until its 
first contact with a wall. In the case of the inactive dimer 
and the pusher motor, once the dimer has reached a wall, 
it remains close to the wall for the remaining time. The 
larger N sphere resides at a fixed distance from the wall. 
This effective attraction to the wall is explained by the 
strong solvent depletion force 45 between the wall and the 
sphere. The smaller C sphere experiences a negligible 
depletion force, which allows the sphere to freely explore 
the wall region that lies within the constraints of the 
dimer bond. The degree of repulsiveness of the wall has 
no discernible effect on the motion of the inactive dimer 
or the pusher motor after being trapped by the wall. The 
puller motor on the other hand behaves very differently 
depending on the strength of the wall potential. For the 
weaker wall repulsive potential, the puller motor is at¬ 
tached to the wall most of the time as well. In contrast 
to the inactive dimer and the pusher motor, however, the 
puller motor occasionally detaches and makes transitions 
to the opposite wall or returns to the original wall. For 
the stronger wall repulsive potential, the puller motor 
is no longer trapped by the walls and frequently makes 
transitions between them. 

The different behaviors of the pusher versus puller 
motor are a consequence of the respective direction of 
propulsion. When the dimer bond is perpendicular to 
the wall, the propulsion force of the pusher motor acts in 
the same direction as the depletion force on the N sphere, 
which reinforces the trapping by the wall. The propul- 




FIG. 2. Trajectories of the perpendicular distances of the N 
(red) and C (blue) dimer spheres to walls, Z , versus time for 
an inactive dimer (top), and puller (middle) and pusher mo¬ 
tors (bottom). The left and right panels illustrate a weaker 
and a stronger repulsiveness of the wall-sphere potential, re¬ 
spectively. The separation between the walls is Lz = 20. 


sion force of the puller motor opposes the depletion force, 
which allows the motor to detach from the wall occasion¬ 
ally or frequently depending on the wall repulsiveness. 

The time- and ensemble-averaged 46 probability den¬ 
sities of N and C sphere positions are shown in Fig. 3 
for the stronger wall potential for an inactive dimer, and 
puller and pusher motors. The insets of Fig. 3 show the 
probability densities plotted as cumulative distributions. 
The N sphere density for the inactive dimer and the 
pusher motor exhibits a pronounced peak, which is the 
signature of the trapping. Since the time until first con¬ 
tact with a wall is shorter for the pusher motor than for 
the inactive dimer, the peak magnitude of the N sphere 
density for the pusher motor is slightly larger and ac¬ 
counts for close to 100% of the corresponding cumulative 
density. The N sphere density for the puller motor shows 
a richer structure with damped oscillations that originate 
from the structural ordering of the atomic-like solvent at 
the wall. The C sphere density has damped oscillations 
for all three cases. Far from the wall, the C sphere den¬ 
sity converges to a uniform density for the puller motor 
and decays to zero for the inactive dimer and the pusher 
motor, identified earlier as the confinement by the dimer 
bond. 
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FIG. 3. Probability densities /(Cs) of the distance Cs between 
each dimer sphere S and the nearest wall for an inactive dimer 
(top), and puller (middle) and pusher motors (bottom), for 
a wall separation Lz = 20. The insets show the correspond¬ 
ing cumulative probability density T(Cs). The dotted lines 
indicate uniform density distributions. 


The structural ordering of the solvent also affects the 
average orientation of the dimer relative to the walls. We 
denote with z the bond vector pointing along the propul¬ 
sion direction from the centers of the N to C spheres 
and separate the contributions perpendicular and paral¬ 
lel to the wall, z± and z M , respectively. The probability 
density of z M assumes a uniform distribution due to the 
periodic boundary conditions in the X and Y directions. 
The probability density of z± is strongly nonuniform, as 
seen in Fig. 4, which plots the joint probability density 
f(z LjCn) °f the perpendicular component of the dimer 
bond vector, z ±, and the distance of the N sphere from 
the closest wall, Cn 5 for an inactive dimer, and puller and 
pusher motors. These probability densities were com¬ 
puted by averaging over slabs [Cat — ,Cn + of 
width AC = 0.4 parallel to the walls. 


FIG. 4. The joint probability density /(«_l,Cn) of the per¬ 
pendicular component of the dimer bond vector, z±, and the 
distance of the N sphere from the closest wall, Cn, for an in¬ 
active dimer (top), and puller (middle) and pusher motors 
(bottom), for a wall separation Lz = 20. The color coding of 
the probability densities is indicated in the inset. The dimer 
configurations associated with the peak maxima are sketched 
in the figures. 


The orientation ranges from z± = —d, which corre¬ 
sponds to the dimer being perpendicular to the wall with 
the C sphere oriented towards the wall, to z± = d , which 
corresponds to the opposite orientation. The densities for 
all three cases show damped oscillations as a result of sol¬ 
vent structural ordering close to the walls, with a period 
equal to the diameter of a solvent molecule. Each peak 
corresponds to a likely configuration of the dimer that 
is described completely by the dimer bond orientation 
and the position of the N sphere. The joint probabil¬ 
ity density f(z±, Cn) plots allow one to assign the motor 
configurations that contribute most significantly to the 
peaks of the probability density. These configurations 
are indicated in the figure. 
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FIG. 5. Mean propulsion velocity, (14), of a confined puller 
motor as a function of the distance of its center of mass to 
the closest wall, £, for e w s = 2. The velocities are averaged 
over slabs [£ — of width = 2.5 parallel to the 

walls. The dotted line shows the propulsion velocity of an 
unconfined puller motor. 

For the inactive dimer and the pusher motor, the ori¬ 
entation densities resemble the C sphere position densi¬ 
ties shown in Fig. 3, which, again, reflects the trapping 
of the N sphere by the wall. For the inactive dimer, 
the density converges to a uniform distribution for close- 
to-perpendicular orientations, while for the pusher mo¬ 
tor, the density decays to near zero. This decay follows 
from the propulsion of the pusher motor parallel to the 
wall, which causes the C sphere to be pulled along by the 
trapped N sphere and thereby closer to the wall. Com¬ 
pared to the inactive dimer and the pusher motor, the 
puller motor exhibits a larger number of likely configu¬ 
rations, since the N sphere is not trapped by the wall. 

The effects of confinement on suspensions of hydro dy¬ 
namically interacting pusher and puller swimmers were 
investigated by Hernandez-Ortiz et al. 25 The swimmers 
were modeled as dimers comprising beads linked by a stiff 
spring, subject to a force that gives rise to propulsion. 
Considering only very low dimer densities, the swimmer 
concentration profile as a function of the distance from 
the walls has peaks near the walls. As one moves farther 
from the walls, the concentration falls to close-to-zero val¬ 
ues for pushers and to non-zero values for pullers. These 
findings are consistent with our puller and pusher single 
sphere-dimer results. However, we note that our sim¬ 
ulations include a variety of other effects including the 
self-generated chemical gradients, solvent structure and 
depletion forces, fluctuations as well as hydrodynamics. 


IV. MOTOR DYNAMICS IN CONFINED GEOMETRY 

The catalytic reaction at the C sphere is responsible for 
a propulsion force along the dimer bond, which, for our 
choice of interaction parameters, results in a mean motor 
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FIG. 6. Parallel orientational correlation function of a con¬ 
fined inactive dimer, and confined puller and pusher motors, 
for Lz = 20 and e w s = 2. The solid lines show fits to an ex¬ 
ponential decay. The inset shows the fitted relaxation times 
r" as a function of Lz- The dotted lines indicate r r for the 
respective bulk case. 


velocity, (14) A 0. The propulsion velocity of a puller 
motor is shown in Fig. 5 as a function of the distance 
from the wall for various wall separations, Lz- Close to 
the wall dimer propulsion is suppressed leading to smaller 
mean velocities of (14) ~ 0.023. Far from the wall the 
velocity converges to a limit value (V z ) ~ 0.029, which 
corresponds to the propulsion velocity of an unconfined 
puller motor. A dimer confined in a narrower slab with 
Lz < 20 has a slightly smaller velocity close to the wall 
than a dimer confined in a wider slab, but overall the 
velocity profile is roughly independent of Lz- 

The dimer is subject to strong thermal fluctuations 
from the solvent that affect its translational and rota¬ 
tional motion. In particular, rotational Brownian mo¬ 
tion will reorient the dimer motor so that its ballistic 
motion will only be manifested for times less than the 
average reorientation time. (In fact, for these Angstrom- 
size dimer motors, the ballistic regime is dominated by 
thermal inertial effects. 32 ) For times longer than this, the 
behavior will be diffusive but with an enhanced diffusion 
coefficient. As discussed earlier, these reorientation ef¬ 
fects are especially strong for Angstrom-size dimers and 
the crossover from ballistic to diffusive motion occurs on 
the time scale of picoseconds. For a dimer motor con¬ 
fined by two parallel walls, rotational motion is no longer 
isotropic and it is interesting to examine how confinement 
influences this property. 

If, as earlier, we denote by z the unit vector along the 
dimer bond, we may resolve this quantity into its com¬ 
ponents perpendicular and parallel to the wall, z± and 
z„, respectively. Here we focus on the parallel component 
since the corresponding correlation function may be com¬ 
pared directly with that of an unconfined dimer. The par¬ 
allel orientational correlation function, C„(t) = (z M (f)-z M ), 
is plotted in Fig. 6 versus time for an inactive dimer, and 
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puller and pusher motors. The parallel orientational cor¬ 
relation follows an exponential decay, C„(t) oc e _4/,T r, and 
the fitted relaxation times t" are shown in the inset of 
Fig. 6 as a function of the wall separation Lz■ For the 
inactive dimer and the puller motor, the confined and 
unconfined dimers exhibit roughly the same reorienta¬ 
tion time. For the pusher motor, however, the relaxation 
times of the confined dimer are almost twice as large as 
for the unconfined dimer. The confinement between the 
two walls has a negligible effect on the parallel orienta¬ 
tional correlation in the case of the inactive dimer and 
the puller motor, but a strong effect in the case of the 
pusher motor. We see this effect for both values of e w s- 

As discussed above, one of the major, easily obtainable 
signatures of propulsion for very small motors is the ex¬ 
istence of enhanced diffusion. Confinement will influence 
this transport property. We have shown earlier 32 that 
the MSD of an unconfined Angstrom-scale sphere-dimer 
motor exhibits an enhanced diffusive regime and is ap¬ 
proximated well by the equation (d is now the number of 
dimensions) 

A L 2 0 1 ) = 2dD m t - 2{V z ) 2 t 2 (l - e” 4 ^) (5) 

Here r v is the decay time of the velocity fluctuations 
as determined from the diffusion coefficient of an inac¬ 
tive dimer, D a = (knT/M m ) r v . The MSD reduces to 
A L 2 (t) « (dk^T/Mm + (V z ) 2 )t 2 in the ballistic regime, 
t <C t v , and to A L 2 {t) ~ 2d(D 0 + g(I4) 2 r r )f = 2dD m t 
in the diffusive regime, t r r . 32 

The MSDs of a confined inactive dimer, and confined 
puller and pusher motors are shown in Fig. 7. To 
account for the confined geometry, the MSDs are sep¬ 
arated into components perpendicular and parallel to 
the walls, AL 2 (t) = ((R±(t) — R±) 2 } and |A L 2 (t) = 
!(|R m (£) — R,,| 2 ), where R (t) is the center of mass of the 
dimer. For the inactive dimer and the puller motor, the 
parallel components show negligible dependence on Lz 
and overlap with the one-dimensional MSD of the respec¬ 
tive bulk case. For the pusher motor, the parallel compo¬ 
nents also show negligible dependence on L z but deviate 
significantly from the MSD of the bulk case. In con¬ 
trast to the parallel components, the perpendicular com¬ 
ponents show L^-dependent behavior. For the inactive 
dimer, the perpendicular MSDs reach a diffusive regime 
for intermediate times for all but the smallest Lz = 10, 
and the diffusion coefficient decreases with decreasing 
Lz- For the puller motor, the perpendicular MSDs show 
the same qualitative behavior as the bulk MSD in the 
intermediate time regime; in the long-time limit, how¬ 
ever, where the bulk MSD reaches the enhanced diffu¬ 
sive regime, the confined MSDs instead converge to Lz- 
dependent limits, which are in good agreement with the 
theoretical limit gA L 2 Z of confined diffusion for a one¬ 
dimensional random walk. 4 ' For the pusher motor, the 
perpendicular MSDs show superdiffusive behavior in the 
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FIG. 7. One-dimensional mean square displacements parallel 
and perpendicular to walls depending on wall separation Lz 
for an inactive dimer (top), and puller (middle) and pusher 
motors (bottom). In each panel, the ensemble of overlap¬ 
ping curves at the top corresponds to the parallel MSDs, 
and the diverging curves below correspond to the perpen¬ 
dicular MSDs. The solid black line corresponds to the one¬ 
dimensional MSD for an unconfined dimer. The dotted lines 
show the ballistic and diffusive regimes of the one-dimensional 
MSD for an unconfined dimer. The dashed, horizontal line in¬ 
dicates the theoretical limit of confined diffusion for Lz = 50. 

intermediate time regime; compared to the puller motor, 
however, the superdiffusive regime is strongly suppressed, 
and in the long-time limit, the MSDs become subdiffu- 
sive, to the point of a plateau for the smallest Lz = 10 
that is two orders of magnitude below |A L 2 Z . The differ¬ 
ences in the L^-dependence of the intermediate regimes 
of the perpendicular MSDs for the three cases follow from 
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TABLE I. Dynamical properties for an inactive dimer, and 
puller and pusher motors, for various wall separations Lz and 
e w s = 2, and the bulk case: Mean propulsion velocity, ( V z ), 
relaxation time of the parallel orientational correlation, t) 1 , 
diffusion constant of the parallel mean-square displacement, 
D'm, and deviation of the theoretical estimate for the diffusion 
constant, D'^ th = Do + 4(I4) 2 t", where d is the effective 
number of dimensions of the system that minimizes th . 


FIG. 8. Mean first-passage time (top) and number of passages 
(bottom) of a puller motor between opposite wall regions as a 
function oi Lz, for e w s = 2. The center of mass of the dimer 
is considered to be in a wall region when the distance from 
the wall £ < 2.5. The right panels show fits to power laws, 
(t p ) oc Lz and N p oc L~^. 


the different structural ordering of the dimer. As noted 
before, the puller motor frequently makes transitions be¬ 
tween the walls for e w g = 2 , which results in long peri¬ 
ods away from the walls with superdiffusive motion, and 
short periods close to the walls with suppressed motion 
perpendicular to the walls. The inactive dimer and the 
pusher motor, on the other hand, are trapped for very 
long periods (as long as our simulation times) after their 
first contact with a wall; their perpendicular MSDs re¬ 
flect the averaging over the initial period of diffusive or 
superdiffusive motion, respectively, and the subsequent 
period of suppressed motion perpendicular to the walls. 
The time until first contact grows with Lz and is larger 
for the inactive dimer than the motor dimers due to the 
absence versus presence of propulsion, which explains the 
more strongly suppressed perpendicular MSDs for the 
pusher motor compared to the inactive dimer. 

The diffusion constants extracted from the parallel 
MSDs of confined dimers and the MSD of unconfined 
dimers are listed in Table I; along with the motor ve¬ 
locities, the reorientation times parallel to the walls, and 
theoretical estimates for the diffusion constant from these 
quantities. For the inactive dimer and the puller motor, 
the diffusion constants D ^ of the confined and uncon- 
fined dimers are approximately equal; for the pusher mo¬ 
tor, the diffusion constant of the confined dimer is twice 
that of the unconfined dimer. The effective number of 
dimensions d in the theoretical estimate for the diffusion 
constant, D }„ = D 0 + ^(V z ) 2 t", is chosen as 2 for the 


pusher motor, which is trapped by a wall most of the 
time, and 2 or 3, depending on Lz, for the puller mo¬ 
tor, which for small Lz spends most of the time in the 
wall regions. For the puller motor, the theoretical values 
of D'^ are in good agreement with the measured values. 
For the pusher motor, the theoretical values consistently 
underestimate the measured values, which, however, is 
also the case for the unconfined pusher motor; hence, 
this deviation is not specific to the confinement. Since 
the reorientation times for the pusher motor deviate sub¬ 
stantially from those of the inactive dimer, this indicates 
that reorientation also has an active component. The 
simple estimate for diffusion enhancement assumes sim¬ 
ple inactive reorientation and this is likely the origin of 
the quantitative failure of the simple theoretical estimate 
for the pusher motor. Nevertheless, these considerations 
suggest that the increased r ) 1 of the confined pusher mo¬ 
tor is responsible for the increased D'^ compared to the 
bulk case. 

For the puller motor that makes frequent transitions 
between the walls, the mean first-passage time and the 
number of passages are plotted as a function of Lz in 
Fig. 8 . The first-passage time and the number of passages 
are found to follow power laws, (r p ) oc L’g and N p oc 
L~p , respectively, with exponent 7 ss 1 . 6 , which indicates 
superdiffusive transitions between the walls . 48 
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V. CONCLUSION 

Molecular motors effect active transport in the cell and 
carry out a variety of other cellular functions. It is cer¬ 
tainly interesting to imagine that, in the future, synthetic 
molecular-scale motors and machines could also be used 
for biological functions at the cellular level. The fact 
that the dynamical properties of small catalysts and sin¬ 
gle enzyme molecules are modified by chemical activity 
supports such a notion. 27-29 While the potential applica¬ 
tions of micron-scale motors are being actively explored 
and many of their properties can be modeled through 
continuum theories, the dynamics of Angstrom-size syn¬ 
thetic motors involve some additional considerations that 
can only be captured through full molecular dynamics 
simulations. They are dominated by fluctuations, molec¬ 
ular reorientation often occurs on picosecond time scales 
limiting the regime where ballistic motion dominates, sol¬ 
vent sizes are often comparable to those of the motor, so 
structural effects come into play and solvent depletion 
forces are very strong. These features, combined with 
the fact that the utility of a continuum description of 
the solvent velocity flow fields must be confirmed, make 
these active objects challenging to study. 

On the Angstrom scale, specific characteristics of the 
molecular nature of the motor and its environment be¬ 
come important. This implies that studies of real sys¬ 
tems on these microscopic scales will have to account 
for these molecular details. The model system studied 
in this paper provides insight into some of the qualita¬ 
tive characteristics that such systems may exhibit. In 
particular, the results show how the distinctive features 
of Angstrom-scale motors, in combination with confine¬ 
ment, change motor spatial structure and dynamics. We 
have seen that self-propulsion is able to counteract sol¬ 
vent depletion forces in some circumstances so that the 
structural ordering of inactive dimers between the walls 
is very different from that of active dimer motors. The 
strong anisotropy and spatial dependence of dynamical 
properties including motor velocity, orientation, and dif¬ 
fusion have their origin in the confining geometry. Such 
information should prove useful for applications involv¬ 
ing very small motors in cellular or microfluidic environ¬ 
ments. 
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